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QQ ■ Abstract 

(N : 

It has been found in in vitro experiments that cytoskeletal filaments 
driven by molecular motors show finite diffusion in sliding motion even 
in the long filament limit [Y. Imafuku et al, Biophys. J. 70 (1996) 
O*! 878-886; N. Noda et al, Biophys. 1 (2005) 45-53]. This anomalous 

Q I fluctuation can be an evidence for cooperativity among the motors in 

■ action because fluctuation should be averaged out for a long filament 
if the action of each motor is independent. In order to understand the 
nature of the fluctuation in molecular motors, we perform numerical 
simulations and analyse velocity correlation in three existing models 

^ ■ that are known to show some kind of cooperativity and/or large dif- 

'sj" ! fusion coefficient, i.e. Sekimoto-Tawada model [K. Sekimoto and K. 

^ Tawada, Phys. Rev. Lett. 75 (1995) 180], Prost model [J. Prost et 

■ al, Phys. Rev. Lett. 72 (1994) 2652], and Duke model [T. Duke, 
Proc. Natl. Acad. Sci. USA, 96 (1999) 2770]. It is shown that 
Prost model and Duke model do not give a finite diffusion in the long 
filament limit in spite of collective action of motors. On the other 
hand, Sekimoto-Tawada model has been shown to give the diffusion 
coefficient that is independent of filament length, but it comes from 

^ ' the long time correlation whose time scale is proportional to filament 

. length, and our simulations show that such a long correlation time 

conflicts with the experimental time scales. We conclude that none 
of the three models do not represent experimental findings. In order 
to explain the observed anomalous diffusion, we have to seek for the 
mechanism that should allow both the amplitude and the time scale 
of the velocity correlation to be independent of the filament length. 

keywords: molecular motor, sliding motion, anomalous fluctuation, coop- 
erativity, model simulation 
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1 Introduction 



One of the outstanding problems in molecular processes in living systems has 
been how they achieve reliable action under the influence of overwhelming 
thermal and/or statistical fluctuations. In the case of muscle contraction, 
Huxley[T] has already noticed, in his original work, that collective action of 
many motors produces smooth sliding motion even though action of individ- 
ual motor is highly stochastic. 

Recently, Imafuku and co-workers have done series of experiments on fil- 
ament motion driven by many molecular motors; Focusing on the fluctuation 
rather than the average motion, they have revealed an intriguing aspect of 
co-operativity in the collective action^, El HI E]. 

They performed the in vitro motility assay on an unloaded filament of 
length L, measured the displacement X{t) over the time interval between to 
and to + t, and evaluated the diffusion coefficient D in the sliding motion 
defined by 



where (■ ■ ■) denotes the average over initial time to and samples. If each motor 
exerting the force to the filament is statistically independent, the fiuctuation 
is averaged out as the filament length L becomes longer, namely, as the 
number of motors that interact with the filament becomes larger. It has 
been shown that D decreases in proportion to 1/L for random action of 
independent motors[6j. In the experiments, however, they have found that D 
is not proportional to 1/L but converges to a constant value for large L. This 
means that the motors are not interacting with the filament independently 
but their actions are correlated with each other. 

In order to understand these results, Sekimoto and Tawada have analyzed 
the motion of a cytoskeletal filament driven by protein motors with random 
orientation, and demonstrated that the diffusion coefficient D of its motion 
is independent of the filament length due to the randomness quenched in the 
motor orientation [7j. Later, however, Noda et a/ have found similar behavior 
of D even for the case where Sekimoto and Tawada model is not applicable, 
i.e. the case where the myosins are not random but aligned|5]. 

No other model has been known to show the constant D so far, and origin 
of the observed fiuctuations has not been understood yet. 

In this paper, we study fiuctuations of cytoskeletal filament motion in de- 
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tail for three existing models: Sekimoto and Tawada model[7|, Prost model[Sl 
El [ini [H], and Duke model[T2]; These are known to produce some kind of 
co-operativity[T3]. As a tool to analyze dynamics, we examine the velocity- 
correlation function of the filament sliding motion obtained by numerical 
simulations. Nature of dynamics shows up in detailed feature of the veloc- 
ity correlation function, and the diffusion coefficient can be derived from its 
integration. 

In Sec. 2, some of the basic formulas are introduced in connection with 
the diffusion coefficient and the velocity correlation. Detailed description and 
results for each model are presented with discussions for Sekimoto- Tawada 
model in Sec. 3, for Prost model in Sec. 4, and for Duke model in Sec. 5. Con- 
cluding remarks are given in Sec. 6. 



2 Diffusion Coefficient and Velocity Correla- 
tion 

Imafuku et a/. [21 [3l H] measured the variance of positional fluctuation of the 
filament defined as 

F^it) ^ (iXit) - {Xm"), (2) 



where X{t) is the displacement over the time interval of length t. The average 
displacement {X{t)) is linear in time, and the mean velocity V of the filament 
is determined by 

t/=limm (3) 

They found that F^{t) increases linearly in time. This is an ordinary 
diffusion process, and is characterized by the diffusion coefficient D defined 
by 

D = hm (4) 

In the actual experiments, it has been found that F^{t) behaves as 

F^{t)^2Dt + a, (5) 

with a constant a, which mainly comes from the finite spatial resolution in 
experiments [2J. 
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The diffusion around the average motion comes from the velocity fluctu- 
ations in the sliding motion, and F^{t) can be expressed as 



F^{t) = 2t 




) 



(6) 



in terms of the velocity correlation 



V){v{to + t)-V)), 



(7) 



where v{t) is the velocity at time t and V is the average velocity. The 
derivation of eq.(l6]) is given in Appendix. 

Note that the relation ([6]) is modified when the data is only available at 
discrete times by the step r as in the experiments or our Monte Carlo (MC) 
simulations: In this case, defining of velocity at time tj = rj as v{tj) = 
[X(tj + t) — X{tj)]/T, the diffusion coefficient D is given by 



in the long-time limit (see Appendix). 

The velocity correlation function C„(t) goes to zero when t becomes large 
enough compared to any relevant correlation time, therefore, the second term 
contribution in eq.([6]) becomes negligible in the large t limit, consequently, 
F^(t) increases linearly in time and D in eq.(jl]) can be expressed by 



namely, the diffusion coefficient is given by the integral of the velocity corre- 
lation function. 

Actual measurements are always based on finite time observations, thus 
we define the finite time diffusion coefficient -Dft(t) as the slope of F^{t) at t, 
then we can show 



namely, only the correlation shorter than t contributes to -Dft(t). 

Experimentally determined diffusion coefficient based on the measure- 
ment of variance for the time interval t corresponds to -Dft(t). This should 
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give a good approximation for the diffusion coefficient ([2]), if t is large enougli 
compared to the correlation times. We should be careful about the effect of 
the measurement time t when we interpret the results, because we do not 
know the length of correlation time for the system in advance. Imafuku and 
coworkers estimated the diffusion coefficient D from experimental data by the 
slope of F^{t) at around t ^ 2s for microtubules driven by kinesin[3] t ^ 0.5 s 
for microtubules driven by dynein|3], and t ~ 0.4 s for actin filaments driven 
by myosin [5]. 

In the following, we analyze the diffusion constant D along with the 
velocity correlation. 

3 Sekimoto and Tawada Model 

Sekimoto and Tawada have proposed a simple model to explain dynamical 
fluctuations in the motion of a cytoskeletal filament driven by protein mo- 
tors fixed on a substrate surface|7j; The motor proteins are assumed to be 
aligned at regular intervals q but in random orientation, and a filament of 
the length L slides over the motors in one direction in the presence of ATP. 
The sliding motion is generated by the motors, which are assumed to attach 
to the filament with the rate constant fcb, make a conformational change or 
"power stroke" to generate force, and detach from the filament with the rate 
constant fcub- The motors interact with the filament independently. If the 
linear force law between the filament and the motors is employed, the slid- 
ing distance generated by a power stroke of the width Oj by the i'th motor 
would be tti/Nb, where Nb is the number of the motors that attach to the 
filament at the time of the stroke. The factor Nj^-^ comes from the fact that 
the motors bound to the filament at the time resist the sliding motion. The 
value of stroke width Oj for the i'th motor is always the same, but different 
at random from that of other motors because of its random orientation. 

Sekimoto and Tawada have demonstrated that the model shows L inde- 
pendent diffusion constant. 

3.1 Simulation method 

We simulate Sekimoto-Tawada model by the following MC procedure with 
the time r for one MC step: (a) Pick a motor, say i, at random out of 
the motors under the filament, (b) If the i'th motor is not attached to the 
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Figure 1: Sliding motion of Sekimoto-Tawada model, (a) Typical trajectories 
of filaments with the length L =3(bottom), 6, and 12 (top) jim. The average 
velocity V — 0.25/xm/s is independent of the filament length L. (b) Velocity 
correlation €^{1) for filaments with L =3, 6, and 12 /im. The inset shows 
short time behaviors (0 < t < 0.6 s). (c) Variance as a function of time 
interval for L =3(top), 6, and 12(bottom) /xm. (d) Length dependence of the 
diffusion coefficient D measured from the slope for L/V <t< 2L/V. 

filament, attach it with the probability /cbT, and advance the filament by 
tti/Nb- If the i'th motor is already attached to the filament, detach it with 
the probability k^hT without any motion of the filament. 

One MC step of r corresponds to repeating this procedure times, where 
N is the number of motors that arc capable of attaching to the filament: 
N — L/q. The results do not depend on r as long as r/cb <S 1 and rfcub ^ 1- 

We adopt the rate constant of the motor attachment /cb = 6.3 s~^, that of 
the detachment fcub =14.7 s^^, and the distance between motors q = 42.9 nm. 
The width of the power stroke is chosen out of random numbers with the 
uniform distribution ranging between and 34 nm with the average 17 nm. 
The time step r is taken to be 0.01 s. 
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3.2 Simulation results 

Figure [U^a) shows typical trajectories of filaments with L =3, 6, and 12 /xm. 
The average velocity V is 0.25/im/s, and is independent of L. 

Figure [T](b) shows the velocity correlation Cv{t) for rather long time (0 < 
t < 50 s), and the inset shows the short time behaviors for < t < 0.6 s. 
We find following characteristic behavior in C„(t): (i) At t = 0, the velocity 
correlation Ct,(0) is the variance of the velocity fluctuation, which takes a 
certain positive value, (ii) For t > 0, Cv{t) drops immediately to almost 
zero, but recovers quickly to reach the maximum, (iii) After the short time 
recovery, Cy{t) decays linearly over a long period of time, (iv) The value of 
Ct,(t) scales as to 1/L for both t = and the maximum value described at 
(ii). 

Note that the actual values of instantaneous velocity variance Ct,(0) do 
not have physical meaning because the model assumes instantaneous dis- 
placement by a stroke of molecular motor. 

The immediate drop of C„(t) comes from the fact that the consecutive 
strokes to the filament is random. The short time recovery is due to the 
correlation given by the strokes from the same motors. The time T when the 
maximum correlation is achieved corresponds to the time interval that each 
motor gives successive strokes to the filament, namely T l/k\^-\- l//cub(Ri 
0.23 s). This time does not depend on L. The gradual hnear decay after 
that comes from the loss of correlation due to the fact that new motors come 
into play as the filament moves. The correlation is lost completely when the 
filament proceeds over the distance L and resides on a new set of motors, 
therefore, the time scale for this correlation is L/V . 

The variance of displacement F^{t) is shown in Fig{T|^c) by a solid line for 
L =3, 6, and 12 /im. We see that the slope increases gradually in course of 
time till t ~ L/V . This comes from the slow decay of the velocity correlation 
Cvit) in Fig{T](b) through Eq.([6]), namely, the slope at time t is given by the 
integration of Cv{s) for < s < t. The slope of the variance is proportional 
to 1/L for a fixed short time t <^ L/V , because C^it) oc 1/L for a fixed time 
t < L/V. The slope for long time t > L/V, however, becomes independent 
of L, because the range that C^,(t) is non-zero is proportional to L. 

Fig{T](d) shows the diffusion coefficient D obtained from the slope for long 
time (fitted in the range L/V < t < 2L/V). One can see that D tends to 
a non-zero constant for larger L. As denoted above, however, D should be 
proportional to 1/L if one measures at fixed t <^ L/V in this model. 
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3.3 Discussions 



The key of this model is that each motor gives strokes of its inherent strength, 
simulating the experimental situation where motors are fixed to a substrate 
surface in random orientation and are not aligned to the direction of the 
filament; Each motor always gives the same stroke, although it varies from 
motor to motor. This makes the correlation time in velocity fluctuation 
proportional to L, and results in the L independent diffusion constant. If 
a stroke by a motor changes randomly every time it gives a stroke, the 
correlation time is independent of L, thus the diffusion coefficient D would 
be proportional to 1/L. 

This assumption of quenched randomness in stroke in the Sekimoto- 
Tawada model is introduced in order to simulate the experimental situations 
by Imafuku et a/.|3l H], where microtubules are driven by motors scattered 
randomly on a substrate. Their results that the diffusion coefficient in sliding 
motion of microtubules is independent of L are reproduced by the model. 

The present analysis shows, however, that this does not necessarily justify 
the model because the experimental time scales are not in the range where 
the model gives L independent diffusion coefficient. In these experiments, 
L/V was 0.5 s to 2 s in Ref.[3] and 0.6 s to 2.5 s in Ref.[l], while the time 
t of the measurement was about 2 s in Ref. [3J and 0.5 s in Ref.[3], respec- 
tively. Namely, the latter experiment by Imafuku et a/.[lj were in the range 
t < L/V, i.e., the range where Sekimoto-Tawada model shows the diffusion 
coefficient that decreases with L. Therefore, the results of the model does 
not correspond to the experimental observation by Imafuku et a/.[l]. 

Relevance of quenched randomness to the anomalous diffusion is also 
questionable; Kinesin has been demonstrated to swivel almost freely to ad- 
just itself in any direction to give effective strokes[Tl]. Furthermore, the 
anomalous filament length independent diffusion has been observed even in 
the experiment where myosins are aligned in orientation [5]. 

4 Prost Model 

Prost and coworkers introduced a two-state model of Brownian motor, which 
exhibits cooperative behavior [HI [10] • In their model, motors are attached to 
a backbone with a fixed spacing q, and each motor can take two states: the 
bound state and the unbound state. In the bound state, the motor is strongly 
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Figure 2: Prost Model 

bound to a filament to exert the sliding force on it, and the interaction 
potential between the filament and the motor is periodic with period / along 
the cytoskeletal filament. A simple saw teeth potential is assumed; Each 
tooth of the potential height U consists of a part with a positive slope of 
length a and a part with a negative slope of length h = 1 — a. The unbound 
state represents the weakly bound state, where the interaction potential is 
fiat, and the motor does not exert force on the filament. The transition from 
the bound to the unbound state takes place with the rate cUub only when a 
motor is around the potential minima within a detaching region of size 5x. 
The transition from the unbound to the bound state, on the other hand, 
occurs everywhere with a constant rate uo\^. The inertia of the filament is 
assumed to be negligible, thus the instantaneous velocity of the filament v is 
determined by the balance between the total potential force from the motors 
Fmot and viscous force —NXv, where N is the number of motors and A is the 
viscous resistance coefficient per motor. 

This model has been demonstrated to show a collective motion of smooth 
sliding when a number of motors are attached to a filament [TlJ, in contrast 
with a Brownian motion under a periodic potential for a single motor system. 
If the periodic potential is asymmetric, the filament shows unidirectional mo- 
tion, while bidirectional motion is observed when the potential is symmetric. 

The cooperativity of the motors is evident especially in the symmetric 
case. In this case, there is no reason for the filament to proceed in one 
direction, but once the filament starts moving to one of the directions by 
chance, it tends to keep on moving in the same direction. This can be 
understood from the fact that the transition from the bound to the unbound 
state occurs only when the motors are around the potential minima; The 
motors are kicked out to the unbound state after they go down the potential 
slope exerting the force to the filament, therefore, there are more motors 
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going down the potential than those going up in the bound state. This state 
of directed motion may be regarded as a state with a broken symmetry due 
to the collective effect [H]. 

The direction of motion flips occasionally when the filament length is 
finite, and the interval between the flips gets exponentially long with the 
filament length. A theory analogous to the one for the phase transition in a 
magnetic system has been developed to describe this state. 

4.1 Simulation method 

One Monte Carlo step with the time r consists of the two procedures: (a) 
the transition trials between the bound and the unbound states of motors, 
and (b) the filament advance. In the procedure (a), each of the N motors 
goes through the following transition trial depending upon its state: if the 
motor is in the bound state and it is located in one of the detaching regions 
of size 6x around the potential minima, detach it to the unbound state with 
the probability cUubT- If the motor is in the unbound state, attach it to 
the bound state with the probability tUbT. Otherwise, do nothing. After all 
of the motors go through the above trials, we perform the procedure (b), 
where the total potential force Fmot acting on the motors in the bound state 
is calculated, and then the filament is displaced by tv with the filament 
velocity v = Fmot/iNX). 

The period of potential is taken to be / = 8 nm, and the motor spacing q = 
42.9 nm with small random distribution of the width ibO.Olnm; The small 
distribution is introduced in order to avoid factitious interference between 
the motor spacing and the potential periodicity. The potential is chosen to 
be a/l = 0.2 and U = 20A;bT with fc^T = 4.14pNnm[g]. The transition 
rates are uj~^ = 2 ms and uj^-^ = 25 ms, and the size of the detaching 
region is 6x = 1.6 nm. We also adopt the viscous resistance coefficient 
A = 1.29 X 10-5 kg/s. 

4.2 Simulation results 

Examples of the simulated displacements of motors versus time under no 
load condition are shown for L =3, 6, and 12 /im in FigJSt^a). The average 
velocity is almost independent of the filament length L, but slightly smaller 
for shorter filaments: V =0.239, 0.248, and 0.252 /xm/s for L =3, 6, and 
12 fim, respectively. The detailed structure of trajectories are shown in the 
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Figure 3: Sliding motion of Prost model, (a) Typical trajectories of fil- 
aments with the length L =3(bottom), 6, and 12(top) /xm. The average 
velocity is almost independent of the filament length L, but slightly smaller 
for shorter filaments: V =0.239, 0.248, and 0.253 //m/s for L =3, 6, and 12 
//m, respectively, (b) Velocity correlation C.u(t) for filaments with L =3(open 
square), 6 (open circle), and 12 (filled circle) /im. The inset shows the behav- 
iors for t > s in larger scale, (c) Variance as a function of time interval for 
L =3(top), 6, and 12(bottom) /xm. (d) Length dependence of the diffusion 
coefficient D. The solid line shows a fitting curve a/L-{l + b/L) via constants 
a = 12.72 X 10-^ /xmVs and b = 2.48 /xm. 
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inset, where one can see more fluctuations in a shorter fllament. Figure 
^h) shows the velocity correlation Cy{t), for which we note the following 
characteristics: (i) C„(t) is largest at t = 0, and decays rapidly to reach the 
negative minimum around t ~ 6 ms. (ii) Then, it increases gradually to 
reach the positive maximum at the time t ~ 30 ms, and decays to zero as t 
increases, (iii) The value of C^(t) is proportional to the inverse of the length 
of fllament, or the number of motors. In Fig. El^c), the variance F^{t)^s are 
shown as a function of time. The diffusion coefficient is evaluated from the 
time dependence of the variance by fitting the data for Os < t < 0.8s to 
Eq.Q. The obtained values of the diffusion coefficient D versus the length 
of the filament L are shown in Fig. El^d) with a fitting curve a/L ■ (1 + b/L). 

4.3 Discussions 

The velocity fiuctuation in this model comes from random transition of mo- 
tors between the two states. Consequently, the amplitude of velocity correla- 
tion is proportional to 1/L for large L, refiecting the fact that the transitions 
of each motor is independent. The time scale of the initial decay in C„(t) is 
set by the transition rate tu's, and the correlation time of the positive peak in 
Cv{t) at if: ~ 30 ms is the time that a motor passes the period of the potential 
{l/v = 8 nm / 0.25 nm (ms)^^ ~ 30 ms), thus, both of the time scales are 
independent of L while oc 1/L, therefore, the resulting diffusion constant 
for this model is proportional to 1/L as is seen in Fig|2]^d). 

5 Duke Model 

Duke [12] proposed a model for the myosin mechanochemical cycle, based 
on the "swinging lever arm" hypothesis, and demonstrated that the model 
shows collective behavior of myosins through coupling of ATP hydrolysis with 
conformational change of myosin head. The conformation change is amplified 
by a lever arm to produce a power stroke. The power stroke stretches the 
myosin neck and the actin filament slides as the neck relaxes. 

The ATP hydrolysis cycle is simplified as in FigJH (a); A main stroke 
of the step width d takes place at the transition from A ■ M ■ ADP ■ Pi to 
A ■ M • ADP, then a small stroke of the width S follows at the subsequent 
transition to A ■ M. 

Eliminating a couple of transient states, i.e. A • M and M ■ ATP, this 
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Figure 4: ATP hydrolysis cycles. 



cycle is further simplified to the three-state cycle shown in Figj4] (b), where 
the small step is included in the transition from A • M ■ ADP to M ■ ADP ■ Pi. 

To each myosin in the attached state A ■ M • ADP ■ Pi, Duke assigned the 
displacement x of the neck, which becomes x + d for A ■ M • ADP after a 
stroke. The elastic energy of the myosin neck at the state A ■ M ■ ADP • Pi 
and that at the state A • M • ADP are Kx'^/2 and K{x + respectively, 
with the spring constant of the neck K. 

The binding rate h from M ■ ADP • Pi to A ■ M ■ ADP • Pi with the neck 
displacement x is affected by this elastic energy; The binding rate to the 
state with the neck displacement in the range [x, x + dx] is given by 



K 



2nkBT 



exp 



Kx^ 
2kBT 



dx. 



(11) 



while the unbinding rate fcub for the opposite transition is assumed to be 
constant and independent of x. 

The transition between the pre and post stroke states, i.e. A ■ M ■ ADP ■ Pi 
and A- M-ADP, is fast and the population ratio is assumed to be equilibrated 
as the ratio of the probability distribution 



P (A ■ M ■ ADP) 
P (A • M ■ ADP ■ Pi) 



exp 



AGstr + AEjx) 



(12) 



where AGstr is a change in chemical free energy and 

AE(x) 



-K(x + df - -Kx^ 
2 ^ '2 



is a change in elastic energy. 

Within the transition from A- M-ADP to M-ADP -Pi, series of transitions, 
including ADP release with the small step 5, ATP bind followed by actin 
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detachment, are involved, and its overall rate is represented as 



^ADP = /^ADP exp 



AEsjx) 



(13) 



with 




(14) 



but the transition in the opposite direction is neglected. 

The time scale for the force balance is assumed to be much shorter than 
the time scale for the transition between the states, thus the position of a 
filament is always adjusted to achieve the force balance among the molecular 
motors as soon as one of the motors changes its state. 

Duke investigated how an ensemble of motors generates sliding motion of 
an actin filament against load. The model displays a transition from smooth 
sliding to synchronized stepwise motion as the load becomes high. 

5.1 Simulation method 

Suppose an actin filament of the length L is sliding straight on a substrate 
under load. The myosin molecules are set in array on the substrate with 
spacing q, thus the number of the myosins that can possibly attach is N ~ 



In each MC step of the time r, the following two procedures are iter- 
ated times: (a) transition between attached and detached states and (b) 
equilibration of attached molecules. Detailed procedures are implemented 
as follows; (a) Transition goes by two steps, (a-1) Pick a motor at random. 
(a-2) If the motor is not attached to the filament, then attach it with the 
probability /c°r with its neck displacement x chosen at random according to 
the distribution ffTTl) . If the motor is attached, detach it with the probability 
fcubT or ^ADpT, depending upon it is in A ■ M ■ ADP ■ Pi or A ■ M ■ ADP, 
respectively, (b) Equilibration is achieved by repeating the following steps 
M times, (b-1) Move the filament to the position where the force from the 
motors balances with the load, (b-2) Distribute all the attached myosins 
according to the population ratio (fT2|) . 

The procedures (a) and (b) are repeated N times in one MC step of time 
r; Within (b), the procedure (b-1) and (b-2) are repeated M times to ensure 
that the system position is relaxed to equilibrium [15]. 



L/q. 
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The parameters used in the simulations are foUowings: the spacing be- 
tween myosin molecules is g = 42.9 nm, the rate constant of the motor attach- 
ment = 40 s~^, that of the detachment fcub = 2 s~^ and ^adp = §0 the 
power stroke size d = 11 nm, the stiffness of myosin neck K = 1 pN/nm, the 
thermal energy fc^T = 4.14 pNnm, the free energy gain AGstr = —IQAksT, 
and the subsequent conformational change 6 =0.5 nm. The time step r is 
taken to be 0.001 s and the equilibration iteration M to be 4. 

5.2 Simulation results 

We show two sets of data in FigO the data for the case without load 
(FigISl^al-4)) and with the load F close to stalling, i.e. F =100 pN//xm(FiglS]^bl- 
4)). We examine the behavior near the stalling load because cooperative 
operation of motors has been observed under load[12j, which could result in 
fluctuations peculiar to the model. For each case, we present trajectories, 
velocity correlations, time developments of displacement variance, and the 
filament length dependences of diffusion constant. 

Figures O^al) and (bl) show trajectories for filaments with L =3, 6, and 
12/xm. The average velocity V is almost independent of L: V =2.06, 2.11, 
and 2.14 /im/s without load, and V =1.98, 2.25, and 2.36 xlO^^/im/s with 
load. The insets show detailed structure of trajectories. One can see larger 
fluctuations in shorter filaments. For the cases under the load, there ap- 
pear back and forth movements(the inset of Figin](bl)); Displacement of the 
movement is of the order of a stroke size d, and its time scale is determined 
by the reaction rates. This is a result of cooperative action of motors; With 
more motors attached, they produce force collectively to move forward, but 
detachment of some motors causes backlash, which leads to further detach- 
ment. 

The velocity correlations are shown in Figsl5](a2) and (b2). For both 
of the cases, immediately after large positive instantaneous correlation, the 
correlation is negative for short time and becomes positive around t = 0.01 s. 
The initial negative correlation is due to the backlash caused by detachment 
after a stroke. The time scale for correlation is determined by the reaction 
rates and does not depend on the filament length. Note that the actual value 
of instantaneous correlation is not meaningful in this model as in Sekimoto- 
Tawada model, because filament position is shifted instantaneously during 
the process of equilibration. 

Difference between the two cases is in the amplitude of the correlation. 



15 




Figure 5: Sliding motion of Duke model (al-4) for cases without load and 
(bl-4) with load of 100 pN//im. (al) and (bl) show typical trajectories of 
filaments with the length L =3(bottom), 6, and 12(top) /xm. The average 
velocity V is almost independent of L: V =2.06, 2.11, and 2.14 //m/s without 
load, and V =1.98, 2.25, and 2.36 xlO~^;um/s with load. The insets show 
detailed structure of trajectories. (a2) and (b2) show the velocity correlation 
Cy{t) for filaments with L =3(open square), 6(open circle), and 12(filled 
circle) /xm. (a3) and (b3) show the variances as a function of time interval for 
L =3(top), 6, and 12 (bottom) //m. (a4) and (b4) show the length dependence 
of the diffusion coefficient D determined by the slope of variances. The solid 
lines show the 1/L fit to the data points. 
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First, the amplitude is much larger for the cases with load than those without 
load. The large correlation in the case under load comes from the back and 
forth movement. Secondly, regarding the L dependence, the amplitude is 
larger for shorter filament for the loadless case as in the previous model, while 
it is smaller for shorter filament for the case with load. The larger correlation 
for shorter filament in the loadless case corresponds to the fact that the 
fluctuations in trajectories are larger for shorter filaments as has been seen 
in Figl5](al). On the other hand, the larger correlation for longer filament 
in the case with load corresponds to the fact that the short time back-and- 
forth movement is more regular in longer filament(the inset of Figl5](bl)). 
The fluctuation in the time scale longer than 1 second seems to be larger in 
shorter filament for the case with load, but the correlation in such long time 
scale cannot be seen in Figl5]^b2) because of statistical errors. 

The time dependences of displacement variance are shown in Figsl5](a3) 
and (b3), and the diffusion constant estimated from these are presented in 
Figsl5](a4) and (b4) with the 1/L curves fitted to the data. For both cases, 
the diffusion constant is proportional to 1/L. This is natural for the case 
without load, because the time scales of C^{t) are independent of L while 
the value of C^it) is proportional to 1/L. In the case with load, the large 
value of the correlation does not give a large value of This means that 
the back and forth movement in short time scale does not give net motion, 
and the diffusion comes from the fluctuation in longer time scale, where the 
fluctuation is larger for shorter filament even in the case with load. 

5.3 Discussions 

The diffusion constant decreases as 1/L with the filament length L. In the 
case without load, the situation is simple; the time scale of velocity corre- 
lation is independent of L and its value is proportional to 1/L, which gives 
D oc 1/L. 

In the case with load, the motors on a filament operate collectively, which 
results in back and forth movement in short time scale. This back and forth 
movement becomes more regular and gives larger value of correlation for 
longer filament, but does not results in larger diffusion. The diffusion in 
the sliding motion comes from longer time fluctuations, which is larger for 
shorter filament. 
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6 Concluding Remarks 



Our results of numerical simulations are summarized as follows: (i) Although 
the amplitude of velocity correlation function decreases as 1/L, Sekimoto- 
Tawada model shows the L independent diffusion coefficient D because the 
time scale of velocity correlation is given by L/V ^ thus proportional to L. (ii) 
Prost model shows D cc 1/L because the amplitude of velocity correlation is 
proportional to 1/L and its time scale is independent of L. (iii) Duke model 
also shows D (X 1/L as in Prost model. 

Although the anomalous diffusion in the experiments is observed without 
load, we examined Duke model with load also, because Duke model has 
been denmonstrated to show a synchronous operation under load and we 
expected that filament motion caused by collective operation of motors was 
a possible explanation for large diffusion. However, the collective movement 
by synchronous operation under load in Duke model did not turn out to 
produce large diffusion for a long filament; It gives rather regular back and 
forth movement in short time, but did not give a net motion for longer time. 

In connection with the experimental observation that D is finite in the 
large L limit, the L independent D in Sekimoto-Tawada model is not likely 
to be relevant; In the experimentally determination of diffusion coefficient, 
only the time dependences of the variance shorter than a few seconds are 
used, which means that the velocity correlation shorter than a few seconds is 
relevant, while the correlation time in Sekimoto-Tawada model is of order of 
L/V, which can be much longer than the experimental time scale for a long 
filament. 

In conclusion, by examining the velocity correlation and the diffusion 
coefficient, we have confirmed the existing models cannot account for the 
anomalous fluctuation that the diffusion coefficient remains flnite even in the 
large L region. Our analysis shows that the experimental observations of 
L-independent diffusion coefficient require a mechanism that makes both the 
amplitude and the time scale of velocity correlation function independent of 
L. 
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A Relation between Variance and Velocity 
correlation Function 



We briefly summarize the derivation for some formulas which express the 
variance F^{t) in terms of the velocity correlation function C„(t) following 
Ref. [TB] for both the continuous and discrete time data. 

A.l Continuous time expression 

The position of the filament X{t) at time t is given by 

X(t) =X(0) + f'v{t)dt, (15) 
Jo 

and the average position becomes 

(X(t)) = X(0) + /Vdt. (16) 
Jo 

Thus, the variance ([2]) is given by 

Fm = fdt'fdt"{{v{t')~V){v{t")-V)) 
Jo Jo 

= 2 I'dt' r dt"C,{t' -t"). (17) 
Jo Jo 

Here, we assumed that the velocity correlation Cy{t) = ((f (to) ~ ^)('^('^o + 
t) — V)) is the function of t and does not depend on to. 

If we take new variables s = t' — t" and s' = t', the integration (1171) 
becomes 

F^{t) = 2 f'ds' r dsC,{s) 
Jo Jo 







s'=t J \ 




s' r dsc.,{s) 


- / ds's'C^is') 




Jo 


s'=0 ^0 / 



= 2t (l - C,{s)ds, (18) 
where we have performed the partial integration. 
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A. 2 Discrete time expression 

For the data at the discrete time sequence tj = jr (j = 0, 1, 2, ■ ■ ■), we define 
the velocity fluctuation as vj = vj — V. The position at time tn is given by 
Xn = Y17=i '^i'^^ its average is given by (a;„) = Vnr. Thus we have 

/ n n \ 

2 



\i=l j=l 



(n n i—1 n j—1 \ 

Y.^i+Y.Y1 ^i^j + EE ^i^j ) 
i=l i=2j=l j=2i=l I 

n i—1 n 

= 2EE(^^^^.)^' + E(^^5r^ (19) 

i=2 j=l i=l 

Now, we assume that Cy(ti,tj) = (viVj) depends only on |tj — tj\, namely, 
Cvitijtj) = Cv{\ti — tj\). Then, the flrst term in (IT^ becomes 

n i—1 n i—1 

i=2 j=l 1=2 k=l 

n—1 n n—1 

= E E c'.(4) = E(^-^)c.(4) (20) 

k=l i=k+l k=l 

while we have Yl'i=i{vf) = nC{0) for the second term in f|T9l) . Therefore, in 
the large tk limit, we have from Eqs. fllQI) and fl20|) with the deflnition (jlj) that 



k=l 



(21) 



D 

which is eq. (j8]). 
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